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Abstract 

We present a study of Y(4260) properties using the initial-state radiation process e + e~ — > 
IisrX (4260). The y(4260) resonance is reconstructed in the tt + tt~ J/ip decay mode, using data 
collected by the Belle detector at the KEKB e + e~ collider. We find a significant signal with a 
central mass value of (4295 ± lOlg ) MeV/c 2 and a width of (133 ± 26+g 3 ) MeV/c 2 . We find 
T ee ■ #(Y(4260) -► vr+vr- J/tp) = (8.7 ± l.lloJ) eV. These results are preliminary. 

PACS numbers: 14.40. Gx, 13.66.Bc, 13.25.Gv. 
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INTRODUCTION 



The Y(4260) was originally seen by BaBar as a significant enhancement in initial state 
radiation (ISR) data in the 7i + n~J/if) final state, where it was fitted using a Breit-Wigner 
function with mass (4259 ± 8±§) MeV/c 2 and width (88 ± 23^) MeV/c 2 ; they found T ee ■ 
B{Y (4260) -> TT+n-J/tp) = (5.5 ± l-0±g;|) eV/c 2 lj|]. BaBar have also seen evidence for the 
decay B~ -> K~Y(4:260), F(4260) -> tt+tt" J/V> 0]. The F(4260) has since been confirmed 
by CLEO in direct production e + e~ — > F(4260) using energy-scan data in addition to 
observing a clear signal in -k + tt~J/iI), CLEO presents evidence for a signal in the K + K~ J /if) 
and tt°7!-°J/iI) final states. Recently, CLEO have also presented results using ISR data, where 
they find a V(4260) — > tt + it~J/?P signal jiJ]. The F(4260) coincides with a minimum in the 
hadronic cross-section [5| and a minimum in the D*D* cross section 0; based on a fit to 
BES data, the authors of Ref. Q set a lower bound i3(F(4260) -> Tr+vr" J/^) > 0.6%. Such 
a prominent hadronic transition to J/^ is unexpected for a cc state above the Z)(*).DW 
threshold, and this has led to various other models being proposed. The F(4260) has 



been described as a 4-quark state 
hybrid 0, 0, MM 



discussed in Refs 17 



a molecular state (9L lid!. U, 12 ]. and a quark-gluon 
the possibility that it could be a conventional charmonium state is 



18 



In this paper, we present the results of a search for F(4260) — ► 7i + 7i^J/if> production 
in ISR events. This study is based on a data sample with an integrated luminosity of 
497.1 fb _1 at the T(4S') resonance and 56.1 fb -1 collected 60 MeV/c 2 below the resonance. 
The data were collected with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 
on 8 GeV) collider 

KEKB operates with a peak luminosity that exceeds 1.6 x 10 34 cm _2 s _1 . The Belle 
detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cerenkov 
counters (ACC), a barrel-like arrangement of time-of- flight scintillation counters (TOF), and 
an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return located 
outside of the coil is instrumented to detect K\ mesons and to identify muons (KLM). The 



detector is described in detail elsewhere [20|. Two inner detector configurations were used. 
A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first sample of 
155.5 fb" 1 , while a 1.5 cm beampipe, a 4-layer silicon detector and a small-cell inner drift 



chamber were used to record the remaining 397.8 fb" 1 m 



RECONSTRUCTION 



Monte Carlo and control sample 



The ISR process was simulated using Phokhara 22| and events were then ported to 
qq98 [Hjj], and processed through a simulation of the detector in GEANT 24|. Decays 
y(4260) — > 7r + 7r~ J /if) were simulated using the if>(2S) as a model [23j. Events were also 
generated for the if>(2S), which is being used as a control sample. A data/MC efficiency 
correction was performed using results from Belle particle identification studies. 
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Track Selection 

Candidate events are chosen from a standard Belle data skim for hadronic events: the 
skimming conditions are optimised for e + e~ — > T(45) — > BB, not ISR events. The skimming 
conditions require the presence of at least three charged tracks (N& > 3), an event vertex 
with radial (r<p) and z coordinates within 1.5 and 3.5 cm of the origin, respectively, a total 
reconstructed energy in the centre-of-mass system (CMS) greater than 0.2^ (yfs is the 
CMS collision energy), a z component of the net reconstructed CMS momentum less than 
Q.5y/s/c, a total ECL energy between 0.1 y/s and O.Sy/s with at least two energy clusters, 
and require R2, the ratio of second and zeroth Fox-Wolfram moments, to be less than 0.8. 
(Here, the z axis is aligned with the symmetry axis of the detector, and opposite in direction 
to the e + beam. The e~ beam crosses this axis at an angle of 22 mrad, so the e + e~ axis as 
seen from the CMS is not quite aligned with z: see the discussion of cos(6 l ) below.) 

We select charged tracks with P t > 0.05 GeV/c for further analysis. Muons are identified 
based on track penetration depth and the hit pattern in the KLM system; the selection has 
an efficiency of 88% for muons. Electron tracks are identified using a combination of dE/dx 
from the CDC, ACC information, E/p (E is the energy deposited in the ECL and p is the 
momentum measured by the SVD and the CDC), and track-cluster matching and shower 
shape in the ECL; our selection has an efficiency of 92% for electrons. To recover radiated 
photons from electron candidates, the four-momenta of the closest photon within 0.05 radians 
of the e + or e~ direction was added to the four-momentum of the track in subsequent 
analysis. Charged pions are identified as charged tracks that fail electron and muon selection 
and additionally pass a kaon veto which is applied using energy loss measurements in the 
CDC, Cerenkov light yields in the ACC, and TOF information. The information from these 
detectors is combined to form a K-ir likelihood ratio, lZ(K/ir) = Ck/(£-k + £k), where 
C-k(£-k) is the likelihood that a pion (kaon) would produce the observed detector response. 
Charged tracks with 1Z(K/ir) > 0.6 are vetoed; this selection has an efficiency of 95% for 
pions, and 13% for kaons. 

Event Reconstruction 

We select events with four charged tracks only, then J /if) candidates are reconstructed 
from J/ if> — > £ + £~, where i is either a positively identified electron or a muon; see Fig. [1] 
top and bottom, where the curves shown are fits using the Crystal Ball function J25[. The 
vertical lines indicate the signal and sideband regions, \M(£ + £~) — m,j/^\ G [0,30] and 
[90, 150] MeV/c 2 respectively. The J/if> candidates are then constrained to a common vertex 
and the nominal J /if) mass to improve the momentum resolution. The perpendicular 
distance between the J/if> vertex and the e + e~ interaction point is required to be less than 
100 /im. We then reconstruct F(4260) candidates by combining the J /if) candidates with 
7r + 7r~ pairs with mass M(n + 7T~) > 0.4 GeV/c 2 . The requirement on M(n + 7t~) suppresses 
the contribution of combinatorial background, including misidentified 7 — > e + e~ conversions. 
The recoil mass squared of the F(4260) candidates, 

M r 2 eco ii= (v^-£y) 2 ~Py 2 , (1) 

is required to satisfy |M 2 ecoil | < 1 (GeV/c 2 ) 2 , in accordance with the F(4260) being produced 
via ISR. 
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FIG. 1: Invariant mass spectrum of J/tp — * e + e~ (top) and J/ip — > (bottom) candidates 

(points with error bars), and the fit described in the text (solid line). The vertical lines indicate 
the signal (red) and sideband (green) regions. 



The ip{2S) is used as a control sample and is reconstructed in the same manner as the 
Y (4260) candidates. Squared recoil mass distributions are shown for the sideband subtracted 
Y(4260) (Fig. [2]) and ij){2S) (Fig. [3]), and compared with signal MC distributions, where the 
signal and sideband regions are defined in Table [I] 

TABLE I: Signal and sideband invariant mass regions, gm^s) = 0.012 is from the wider of the 
Gaussians fitted to the ip(2S) signal. 



Region 


Mass range (GeV/c 2 ) 


F(4260) Signal 


4.2 < M(Tr+ir- J/ip) < 4.4 


Y~(4260) Sidebands 


3.9 < M(vr + vr- J/ip) < 4.0 




5.1 < M(ir+Tr-J/ip) < 5.5 


rp(2S) Signal 


|M(7r+vr- J/V>) - 3.686| < 3a = 0.037 


ip(2S) Sideband 


4cj = 0.049 < \M(ir + Tr~J/ip) - 3.686| < 0.085 = 7a 
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FIG. 2: Recoil mass squared for sideband subtracted y(4260) candidates in data (points) and MC 
(histogram). 




FIG. 3: Recoil mass squared for sideband subtracted ip(2S) candidates in data (points) and MC 
(histogram) . 
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FIT 



We perform a binned maximum likelihood fit to the distribution of invariant mass m = 
M(ir + 7T~ J/ip). The signal function is the product of a Breit-Wigner function, a phase 
space term, and a mass-dependent efficiency correction, added to a second-order polynomial 
representing the background. The Breit-Wigner term is 

1 

(m - M Y ) 2 + ir 2 ' 

the phase space term is 

V /(K + - M 2 +7r _)/(2m)) 2 - m^, 

and the correction for the efficiency as a function of mass, based on an interpolation of 
results from F(4260) MC samples generated with different central mass values, is 

a e ■ (m - M ) + b e , 

where M = 4.3GeV/c 2 ; we fix the parameter M n + n - = 0.5GeV/c 2 . We have neglected 
the effect of interference with other resonances. The slope and intercept were determined, 
respectively, to be a e = (7.4 ± 1.3) (GeV/c 2 ) -1 and b e = (9.31 ± 0.07). The parameters 
allowed to float in the fit were the number of signal (N) events, the mean and width of the 
Breit-Wigner function (My, T), and the background parameters. 




FIG. 4: Invariant mass for y(4260) candidates in data (points) and J /ip mass sidebands (shaded 
histogram), with the fit to data (solid line) and its background component (dotted line). 

The fit to the data is shown in Fig. HI together with the (featureless) J / ip mass sidebands 
(shaded histograms). The signal parameters were found to be M Y = (4295 ±10^°) MeV/c 2 , 
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FIG. 5: Invariant mass for ip(2S) candidates in data (points), with the fit to data (solid line) and 
its background component (dotted line). 



r = (133 ± 26lg 3 ) MeV/c 2 with a yield N = (165 ± 24+£ 3 ) events, where the systematic 
terms will be described in the next section. The significance of this signal was estimated by 
comparing fits with and without a signal term and was found to be Her. 

The ip(2S) control sample was fitted using two Gaussians and a bifurcated Gaussian, 
with a common mean M^s) floating in the fit, and a linear background: see Fig. [51 Using 
the Particle Data Group value of T ee (tp(2S)) jl] we expect our ip(2S) yield to be 4575 ±290, 
where the error is dominated by MC statistics. We observe a ip{2S) yield of 4188 ±67, where 
the error is statistical only. 



ISR photon reconstruction 

Photons are identified as ECL energy clusters that are not associated with a charged 
track and have a minimum energy of 0.060 GeV. If the photon with highest energy E 1ISR 
in an event satisfies E 1ISR + E^+^-j/^ > 10 GeV then the ISR photon is considered to 
be reconstructed. The F(4260) candidates that pass this criteria are shown in Fig. [6J an 
enhancement is seen at M(tt + tt~ J/ip) ~ 4300 MeV/c 2 . 



PROPERTIES 
Dipion Mass 

Figure [7] shows the efficiency-corrected dipion mass distribution for F(4260) candidates 
with the prediction from MC, based on fits to M(tt + 7t~ J/tp) in M(7r + 7r~) bins. The effect 
of the 7t + 7r~ J/ip mass threshold that is introduced by the M(7r + 7r~) binning is taken into 
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FIG. 6: 1"(4260) invariant mass spectrum requiring reconstruction of the ISR photon. 




M(7i + 7i") (GeV/c 2 ) 



FIG. 7: Invariant mass of tt + tt combinations in y(4260) candidates from fitted yields, after 
efficiency correction in data (points) and MC (histogram). 



account by multiplying the signal and background terms by an appropriate threshold func- 
tion in each bin. Both signal and background parameters are fixed from the fit to the full 
sample, and only the normalisation and threshold function terms are allowed to vary. The 
yield (ni) in each bin is corrected using the efficiency (e*), calculated in the same bin using 
F(4260) MC. 
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FIG. 8: Invariant mass of tt + tt combinations in ip{2S) candidates, after efficiency correction in 
data (points) and MC (histogram). 

The ip(2S) dipion mass distribution in Fig. [8] is corrected using a 2nd order polynomial 
fit to the binned efficiency from MC produced at the nominal ip(2S) mass, the errors are 
dominated by the MC error in the binned efficiencies. 



cos(#) distribution 

Figure [9] shows the efficiency-corrected cos(#) distribution for F(4260) candidates with 
the prediction from MC, where 9 is the polar angle of F(4260) candidates relative to the 
direction of the e~ beam, as seen from the CMS. M(tt + tt~ J/ip) fits are performed in cos(#) 
bins, with signal and background parameters fixed from the fit to the full sample, and 
only the normalisation terms allowed to vary. The yield in each bin is corrected using the 
efficiency calculated in the same bin, based on the F(4260) MC as above. Smaller bins were 
used close to cos(#) = ±1, where the efficiency is changing rapidly, and then combined for 
the final plot. The same procedure was followed to produce the cos(6 l ) plot in the ip(2S) 
reference channel, Fig. [TUJ 
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FIG. 9: Efficiency corrected cos(#) distribution for Y(4260) candidates in data (points) and MC 
(histogram). 
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SYSTEM ATICS 



The systematic uncertainties on the mass and width measurement were estimated by 
varying the fitting procedure. Variant fits including signal functions without phase space 
and efficiency corrections, a linear background shape, 10 MeV/c 2 binning, and a (3.8,5.0) 
GeV/c 2 fit range were performed, and the largest deviations were taken as the positive and 
negative systematic errors due to fitting. We then add in quadrature terms from possible 
biases in F(4260) mass reconstruction (based on a fit to MC) and the overall mass scale 
(using the if>(2S) for calibration). The fitting function does not take the variation with mass 
of the ISR cross-section into account. We estimate the resulting error in T ee ■ B (see the 
next section) using the relative change in the second-order QED radiator W(s, x) [26j as the 
7T + 7T~J/ip mass is changed from the fitted mean value by ±1T. (x = 2E~ f /y/s, where E 1 is 
the energy of the ISR photon in the CMS.) A summary of systematic error terms can be 
found in Table [III 



TABLE II: Sources of systematic error. 



Source 


Mass 


Width 


a-B 


r ee -£ 




(MeV/c 2 ) 


(MeV/c 2 ) 


(fb) 


(eV) 


Fitting procedure 


+10 

-2 


+13 
-6 


+i 

-5 


+0.2 
-0.9 


Measurement of mass in y(4260) MC 


±2 








Mass scale from measurement of ip(2S) 


±2 








QED radiator at masses of ±ir 








±0.2 


Total 


+10 

-3 


+13 

-6 


+1 

-5 


+0.3 
-0.9 



Additional fits were conducted on the di-electron and di-muon samples independently 
and found to be consistent with the combined result. 

CROSS SECTION 

Taking the fitted Y (4260) yield 7ij and the efficiency e, in bins of dipion mass, as described 
above, we identify 

— = fdtC ■ a(e + e~ -> 7/5R F(4260); = 10.58)-£(F(4260) -> n+n-J/iP)- 

b(j/i/> -> tr). 

The use of M(ir + ir~) bins reduces the model dependence of the result at the cost of increasing 
the statistical error. Substituting the world average value for B(J/i/j -> £ + £~) @, we find 

a(e + e~ -»■ jisrY; ^/s = 10.58) ■ i3(F(4260) -»■ tt+tt" J/ip) = (48 ± 6+^) fb. 

Using equation (7) from Ref. [26|] we calculate, in the small width approximation, T ee ■ 
£(F(4260) -> Ti+Ti- J/if>) = (8.7 ± l.l±g;|) eV. 

From the signal function of a fit to the full sample (without applying an M(tt + 7t~) cut) we 
find the cross-section at the peak value of 4.295 GeV/c 2 for e + e _ — > F(4260) — ► ti+ti^ J /ip 
to be 43 pb. 
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DISCUSSION 



We have confirmed the y(4260) — > 7t + tt~J/iIj signal with a significance of 11a. The 
dipion mass distribution favours high values of M(tt + 7t~), which is consistent with BaBar's 
result; there is also some evidence for a rise in cross-section near the M(ir + ir~) threshold, 
although statistical errors are large. The cos(6 l ) distribution peaks strongly at values near 
±1, consistent with the interpretation that the events are due to ISR, and when we require 
the reconstruction of the ISR photon we see an enhancement near the central mass value. 

For our signal we find a mass of (4295 ± 10±3°) MeV/c 2 . There have been predictions 
that a hybrid would favour decays to an S-wave plus a P-wave state j27|; the relevant mass 
threshold is M(D°D ( () = 4287MeV/c 2 . We find a mass above this value, although the 
difference is not statistically significant. CLEO's preliminary measurement is (42831 \l ± 
4) MeV/c 2 , coinciding with the threshold and consistent with our value; BaBar's published 
mass is (4259 ± 8±1) MeV/c 2 , below threshold and 2.6a below our value. Our measurement 
of the width, T = (133±26^ 3 ) MeV/c 2 is consistent with those of CLEO (70+^±5) MeV/c 2 
and BaBar (88 ± 23^) MeV/c 2 ; errors are large. Our result for T ee ■ B{Y — > it + it~ J/ip) is 
about 1.6a higher than BaBar's reported value (5.5 ± LOl^f) eV/c 2 . 

SUMMARY 

We have investigated the properties of the F(4260) using the initial-state radiation process 
e + e~ — > 'JisrY (4260). We observe a significant signal described by a Breit-Wigner with 
mass (4295 ± 10t\°) MeV/c 2 and width (133 ± 26^ 3 ) MeV/c 2 , and find T ee ■ £(F (4260) -> 
n + n~J/ip) = (8.7 ± l.llg^eV. These results are preliminary. 
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